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Abstract 

Assessing the significance of apparent total ozone trends is 
equivalent to assessing the standard error of the means. Standard 
errors of time (area) averages depend on the temporal (spatial) 
variability and correlation of the averaged parameter. Trend 
detectability is discussed, both for the present network and for 
satellite measurements, using statistics from daily observations at 
Dobson stations from 40° and 60°N. 


1. Introduction 

For several years much interest has been attached to detection 
of possibly anthropogenic trends in total ozone, either at single 
stations or station groups. Significance of trends or, equivalently, 
the standard errors of point- or area-means, is properly derived 
from knowledge of variances and of data independence, i.e., knowledge 
of temporal and spatial autocorrelations (e.g., Lieth, 1973; Jones, 
1975). In general, authors who report ozone trends (e.g., Angell 
and Korshover, 1973, 1976; Komhyr et al . , 1971, 1973; London and 
Kelley, 1974; Hill et al . , 1975, 1977) use only monthly mean data, 
and are not explicit about how they assessed the standard error of 
the monthly averages, or, where used, of the area averages. The 
purpose of this note is to present estimates of standard errors of 
total ozone time and area means, as derived from ozone's natural 
temporal and spatial variability and autocorrelation in middle 
latitudes determined from daily Dobson data. The use of this in- 
formation in assessing detectability of total ozone changes, at singl 
stations and over areas, will be demonstrated. 



2. Method 


a. Data . Daily total ozone data for each of 26 Dobson stations 
between 40° and 60°N for the period 1957-1972 were obtained from the 
World Data Center for Ozone, Toronto, and were checked for gross 
errors before processing. For our purposes, a trend is defined as a 
change, of time scale at least one year, which is not explained by 
deterministic variations. Trends must be detected against a backdrop 
of non-determini Stic variability, and a primary task is to describe 
this variability. In order to do this, the mean, a trend over the 
entire period-of-record, an average 29-month quasi-biennial oscillation, 
and the first three harmonics of the annual variation were subtracted 
from the data (see Wilcox et al., 1977) and the study proceeded using 
the residuals. These residuals primarily contain a somewhat per- 
sistent ("reddish") synoptic scale variability, but with assumed 
"white" contributions from smaller scale processes and from obser- 
vational error. Any unremoved deterministic periodicity will in- 
crease the correlations, but this effect is thought to be quite small. 
The residuals are undoubtedly more seriously affected by slow 
calibration drifts and by changes in wavelengths used, and there 

needs to be a comprehensive, continuous program to check calibrations 
and observation techniques, as well as to recompute published values 
as necessary. The subtracting of a trend from the data can help 
remove slow calibration drifts, and it is hoped that any remaining 
non-random observation error is relatively small. 

b. Standard error of time averages . For an atmospheric variable 
whose autocorrelation, R, is approximated by a "red-noise" model 6 = 
exp(-bx), where x is lag, Lieth (1973) has shown that 

2 

” -CTtl-exp(-bT)]) 

a 

Here, Oj is the standard error of the mean, a the standard deviation 
of the unaveraged time series, and T the averaging interval. 

To find this ratio for total ozone, correlations at lags from 1 
to 16 days were computed from the residuals for the four seasons 
(winter is December through February). "Zonal mean" correlation 
coefficients were estimated by weighting the station correlation 
coefficients by the square root of the number of observation pairs 
at that station relative to the total pairs of all 26 stations. 

The resulting correlation coefficients for lags 1-7 days are shown 
in Fig. Al. The average number of pairs for any lag at a station was 
465 in winter and slightly larger in the other seasons. For the 
purpose of assessing the significance of the zonal mean autocorrela- 
tions it is estimated that, of the 26 available stations, nine are 
independent (more on this in the next section). Effectively, then. 
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Figure Al. Total ozone temporal autocorrelations, zonal ly 
averaged, from the Dobson stations between 40° and 60°N, 
for lags 1 to 7 days. 



we have about 4200 pairs at each lag, which implies that any "zonal 
mean" correlation coefficient above about 0.04 in absolute value 
is significant at the 99% level. 

The computed coefficients for lags 1-5 days were fitted with 
the simple red-noise model 6 (t) = a^exp{-b t). The coefficients a 
and b , for each season, are shown Tn Tabl5 Al. The "zonal mean" ^ 
temporal standard deviations, both day-to-day, aj, and year-to-year, 
cra> were also obtained by weighting the variances at stations by the 
square root of the number of observations. 

The difference of a^- from unity is mostly due to observation 
error (see e.g., Julian & Thiebaux, 1975), and the value 
2 1 2 2 

=2 (1-a^ )a^ , whose square root is given in Table Al, may be 

interpreted as variance due to observation error. Therefore, 
reduced by this amount is the "true" day-to-day variance. Using 
the values of b.^ from Table Al, equation (Al) yields oj/a = .21, .21, 

.21, and .22 for winter, spring, summer, and fall, respectively 
(90-day means). For example, a single winter season average at a 
typical mid-latitude station will have a standard error (from the 
true season mean) of .21a, or 8.2 D.U. Calculations for other 
averaging times and/or seasons can easily be made. 

It is a well-known result of sampling theory that 

= a^/N (A2) 

where N is the effective sample size. Since the asymptotic value of 
2 2 

ay /a , from equation (Al), is 2/bT-T, the effective sample size N = T/Tq 
where oo 

Tq = / R(T)dx = 2/b^ (A3) 

— oo 

is a characteristic time between effectively independent observations 
(Lieth, 1973). Values of T^ are shown in Table Al. 

c. Standard error of area averages . The preceding method may 
also be applied to the determination of errors in spatial averages and 
equation (Al) again applies, with T now understood to refer to an averaging 
distance. Correlation coefficients of nearly simultaneous residuals 
at different stations were computed, and, in order to include large 
lags while keeping some uniformity in the data set, pairs were 
restricted to those whose orientations were more east-west than 
north-south. Figure A2 shows these correlation coefficients (to save 
space, only those pairs whose longitude separations are 90 degrees or 
less are shown). Similar computations have been made by Fabian (1967) 
for a set of European Dobson stations and by Nastrom (1977) for air- 
craft measurements of ozone concentration near the tropopause. 
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Table Al. Total ozone temporal variability and correlation statistics, zonally averaged 
from Dobson stations between 40° and 60°N. 


‘^d Time 


Season 

Model 

a 

T 

T 

b 

T 

"Day-to-day" 

standard 

deviation 

^d 

"Year-to-year" 

standard 

deviation 

‘^a 

observation 

error 

a 

e 

corrected 

for 

obs. error 
a 

between 

independent 

observations 

To 

Dec-Feb 

.940 

.490 

40 D.U. 

15 D.U. 

9.7 D.U. 

39 D.U. 

4.1 days 

Mar-May 

.856 

.470 

36 

11 

13.1 

34 

4.3 

Jun-Aug 

.812 

.493 

22 

9 

9.2 

20 

4.1 

Sep-Nov 

.880 

.448 

25 

9 

8.4 

23 

4.5 

Annual 

Average 

.872 

.475 

31 

11 

10.1 

29 

4.3 


CORRELATION COEFFICIENT , 


R=l.04 exp(-0.078 \) 


R= 1.05 exp (-0.099 X) 


30 40 50 


R=l.05 exp(-0.l33 X) 


R»l.02 exp(-0.09I X) 


X. LONGITUDE SEPARATION (DEGREES) 


Figure A2. Total ozone spatial autocorrelations for separations 
(lags) 0 to 90 degrees longitude. Only Dobson stations between 
40° and 60°N were used, and the pairs were restricted to those 
whose orientations were more east-west than north-south. 












An average of about 400 observation pairs goes into each 
correlation coefficient in Figure A2. However, the autocorrelation 
within the time series at each station affects the significance of 
the cross-correlation between the two stations. Mitchell's (1963) 
approximation for effective sample size using purely persistent 
series implies that a correlation coefficient must be above about 
0.18 in absolute value to be significant at the 99% level. 

The red-noise model R = a^exp(-bp^x) , where A is longitudinal 
separation (lag), was fitted to various sub-interval averages of 
the measured correlations from the range of lags 0 to 20 degrees. 

In computing the sub-interval averages, the individual correlations 
were weighted by the square root of the number of observation pairs. 
The time difference between observations at both stations of a pair 
was generally less than about 2 hours for these small separations, and 
the associated temporal variability was neglected. (Note that the 
taking of sub-interval averages effectively lowers the 0.18 
significance threshold somewhat in the range 0-20.) The values of a 
being generally greater than 1 indicates that a better model would 
have used A raised to a power slightly greater than 1; however, this 
refinement in the present study does not seem warranted. Inserting 
values of b^^ in equation (Al) yields, for zonal means, a-r/a = .26, 

.23, .20, and .24 for winter, spring, summer, and fall respectively. 
Application of these results will be demonstrated presently. 

An effective length between independent observations, Lg = 2/b^, 
is also given in Table A2. These values were used in subjectively 
estimating, in the previous section, that of the 26 available stations 
only about 9 were independent. Extrapolation would indicate that on 
the order of 100 effectively independent daily values are possible 
from an ideal global network. 


Table A2. Total ozone spatial variability and correlation statistics, 
zonally averaged, from Dobson stations between 40^ and 60° N. 


Season 

Model R = 


Longitude separation of 
independent observations 

^0 

Dec-Feb 

1.04 

.078 

26 deg 

Mar -May 

1.05 

.099 

20 

Jun-Aug 

1.05 

.133 

15 

Sep-Nov 

1.02 

.091 

22 

Annual 

Average 

1.04 

.100 

21 


3. The significance of point- and area-average total ozone trends. 

A common method for determining trends uses monthly or seasonal 
deviations from the long-term normal (e.g., Angell and Korshover, 1973, 
1976; Pittock, 1974; London and Kelley, 1974). Significance of single 
station trends is determined by the size of the standard error of these 
deviations, which can now be determined from equation (Al), using T = 30 
or 90 days, and a from Table Al . If long-term means are required, the 
problem is only slightly different. Here, it seems reasonable to assume 
that little interannual correlation exists in the residuals (e.g.. Hill 
et al., 1975), except possibly for a small amount due to a sunspot 
cycle. Therefore, the standard error olj of a long-term seasonal 
(90-day) mean as an estimate of a climatic mean is given by equation (A3) 
substituted in equation (A2), i.e.. 



(A4) 


where Y is the number of years considered, and the inclusion of a 
accounts for the effect of interannual variability. 

Several authors (e.g., Angell and Korshover, 1973, 1976) compute 
means at groups of stations in order to estimate regional trends. To 
determine the standard error of such group means, one must account 
for both the temporal and spatial correlation. To fix ideas, consider 
the standard error of a group mean for several stations which can 
all be conveniently enclosed in a rectangle whose sides are lengths 
Li and Lo. The standard error of the seasonal mean at each of the 
stations'^is given by equation (Al), and as a first approximation we will 
consider that the mean at every point within the area is known 
within the same standard error. This is an optimistic view, but one 
which becomes more realistic as station density increases. To 
account for the spatial averaging, it is appropriate to apply equation (Al) 
using b. with T = L., and then apply it again using T = Lg. 

This assumes isotropy, which again is probably valid only to a first 
approximation (see, e.g., Buell, 1972; Julian and Thiebaux, 1975). 

Note also that the spatial standard deviation of the (instantaneous) 
field of total ozone is now required. This computation has not been 
carried out; however, oscillations in the residuals are likely due 
predominantly to truly transient eddies, and thus will affect all 
stations in the latitude band more or less equally. This being the 
case, the temporal standard deviation should be a reasonable approxi- 
mation of the spatial standard deviation. 

As an example, the standard error of yearly means for the group 
of North American stations between 40° and 60°Nwill be estimated 
(Churchill, Edmonton, Goose Bay, Caribou, Green Bay, Bismark, Bedford, 

Fort Collins, Boulder, and Toronto). Using the "zonal mean" 
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statistics which are already given in Tabi e A1 (instead of a set 
derived specifically from these stations), oj/a - 0.11. The 
longitudinal extent of this area is about 40 degrees, while the 
latitudinal extent is about 17 degrees, or roughly equivalent to 
26 degrees longitude. Using equation (Al), these values of T yield 
ratios aWa = .61 and .70 respectively, and the result if oj = 

(.11 X .61 X .70)a = .047 (29) =* 1.4 D.U. In other words, a detected 
2oj (i.e., 2.8 D.U.) change from one year to any other year in the 
annual mean total ozone in this region could be judged significant 
at the 95% level of confidence. 

4. Concluding remarks 

It has been shown that standard errors of the mean for time- and 
space-averages are properly determined from time and space correlation 
and variability statistics. Sample statistics have been given, on a 
seasonal and zonal mean basis, for the Dobson stations between 40° 
and 60° N. Further research in this vein should aim at determining 
these statistics for all the specific regions and years where 
Dobson measurements are available. Also, since anisotropy is 
expected, the north-south statistics should be included. 

Hill et al . (1977) , have assessed the detectability of global 
total ozone trends at about 1%, assuming an (independent) 18-station 
network. The present work suggests that there probably exist at 
least that many stations whose monthly means are independent (but 
see Pittock, 1974). However, Hill et al., apparently assume that 
these 18 determine the means (to within the same standard error) at 
all points on the globe. This does not seem likely, as has been 
previously stated implicitly by Kohmyr et al . (1971); Pittock (1974); 
and Angel! and Korshover (1976). However, the example in the 
preceding section points to the likelihood that trends can be 
calculated with high confidence for certain areas of high station 
density, such as the United States and southern Canada, Europe, 
and perhaps India and Japan. Trends over regions not presently 
well -sampled will apparently only be detected by future satellite 
observations of total ozone. Note that the assumption that one 
knows the mean to within the same standard error at all points 
within a region is well-satisfied with satellite observations. 

It is thus of some interest to determine, a pri ori , the space and 
time scales of averaging required to detect a change in ozone with 
a given level of significance from satellite observations. For 
example, let it be necessary to detect a regional 1% change in 
annual mean total ozone (i.e., 3.5 D.U.) at the 99% significance 
level. Such a confidence level requires that 2.6aj < 3.5 D.U., 
or Qj £ 1.3 D.U. Averaging over a year (with b = .475 and o = 29) 
yields Oj = 3.1 D.U. In addition, averaging ov5r a square 3000 km 
on a side would decrease oj/a by another factor of (.60)^ to 
a-j- = 1.1 D.U., thus satisfying the requirement. Such preliminary 
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estimates of trend significance are important in the planning of 
observing programs. However, several years of satellite observations 
will have to be available before year-to-year changes can be 
recognized as either part of some periodicity or not. 
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